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Abstract: We compute the integrand of the full-colour, two-loop, hve-gluon scat¬ 
tering amplitude in pure Yang-Mills theory with all helicities positive, using general¬ 
ized unitarity cuts. Tree-level BCJ relations, satished by amplitudes appearing in the 
cuts, allow us to deduce all the necessary non-planar information for the full-colour 
amplitude from known planar data. We present our result in terms of irreducible 
numerators, with colour factors derived from the multi-peripheral colour decompo¬ 
sition. Finally, the leading soft divergences are checked to reproduce the expected 
infrared behaviour. 
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1 Introduction 

Precision measurements at the Large Hadron Collider provide detailed information 
about the nature of the strong interaction and its role within the Standard Model. 
With new data already arriving from Run II, there is a growing need for higher- 
precision theoretical predictions for a variety of different observables. Recent years 
have seen considerable progress in this respect. Our ability to make predictions in 
perturbative QCD now covers most of the relevant 2 —)■ 2 scattering processes at 
next-to-next-to-leading order (NNLO) precision [1-13], as well as the example of 
inclusive Higgs production at N^LO [14]. Despite this, processes with more than two 
particles in the final state remain beyond the reach of current NNLO methods. 

Next-to-leading order (NLO) corrections to high multiplicity final states are by 
now commonplace in phenomenological studies. Such computations are possible 
thanks to automated techniques, which make use of integrand reduction [15], re¬ 
cursive techniques [16], (generalized) unitarity cuts [17-20] and the known basis of 
scalar integrals. Processes with up to five coloured partons [21, 22] in the final state 
are feasible using on-shell methods that, by only working with the physical degrees 
of freedom, are efficient at controlling the complexity. 
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Multi-loop calculations of 2 —)■ 2 scattering processes [12, 23-32] have been 
quite successful using the more traditional approaches of Feynman diagrams and 
integration-by-parts identities (IBPs) [33, 34], though there are notable exceptions 
using unitarity cutting techniques [35-40]. The main focus for these processes has 
been in the evaluation of the resulting master integrals. At higher multiplicity, rapid 
growth in the complexity of the Feynman diagram representation motivates an al¬ 
ternative approach that makes use of the lessons learned during the automation of 
one-loop computations. 

Two methods have already been explored in this direction. The first of these is 
maximal unitarity [41], which generalises the cutting techniques of Britto, Cachazo 
and Feng [19] and Forde [42] to compute the rational coefficients of the master in¬ 
tegrals, incorporating information from IBPs. Maximal unitarity has been used to 
look at maximal cuts for a variety of high-multiplicity examples in four dimensions 
[43-50]. The second approach extends the integrand reduction program of Ossola, 
Papadopoulos and Pittau [15]. The initial steps in this direction [51, 52] have now 
developed into a deeper understanding using the language of computational alge¬ 
braic geometry [53-57]. The H-dimensional extension of this method has also been 
understood and applied in the context of the planar two-loop, five-gluon amplitude 
in QCD with all helicities positive (all plus) [58]. 

In the context of supersymmetric theories, computational methods based on an 
analysis of unitarity cuts have enabled a large number of high-loop computations. 
Other methods have also been developed, mainly in the context of these simplified 
theories. For example, the colour-kinematics duality of Bern, Carrasco and Johans¬ 
son (BCJ) [59, 60] has been successfully exploited to hnd the complete colour-dressed 
four-loop, four-gluon amplitude in A/" = 4 supersymmetric-Yang-Mills (sYM) [61]. 
The two-loop, hve-gluon amplitude, computed in A/" = 4 sYM in ref. [62], has since 
been extended to the non-planar sector and cast into a complete set of numerators 
satisfying colour-kinematics duality [63]. The integrands in planar A/" = 4 sYM are 
known to all loop orders [64, 65] , and recent studies indicate that this simplicity may 
extend to the non-planar sector [66, 67]. 

The observation that A/" = 4 sYM theory and the all-plus sector of QCD^ are 
related by a dimension-shifting relation [68] suggests that the all-plus amplitude at 
two loops could be a useful testing ground for new techniques. Indeed, at two loops 
the planar sector of the all-plus amplitude was observed to be related to the A/" = 4 
amplitude at the integrand level in a similar pattern to the one-loop story, although 
additional corrections to the A/" = 4 sector appeared in the form of one-loop squared 
(or butterfly) topologies [58]. This fact prompts the question as to how much the 
techniques applied in supersymmetric cases may help to simplify QCD applications. 

In this article, we complete the computation of the two-loop, five-gluon, all- 

^The all-plus sector is equivalent to self-dual Yang-Mills at one loop. 
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plus helicity amplitude including the non-planar sector. In order to deal with the 
increase in complexity of the full colour amplitude, we introduce a method to hud 
compact colour decompositions that make full use of the underlying Kleiss-Kuijf 
(KK) relations [69] in a similar way to the previous treatment at tree level and one 
loop by Del Duca, Dixon and Maltoni (DDM) [70, 71]. We then further exploit 
the on-shell construction of the irreducible numerators to show that all of the non- 
planar cuts can be obtained from the planar cuts. This is reminiscent of the colour- 
kinematics duality, and indeed we employ the BCJ relations [59] at tree level to relate 
planar and non-planar cuts. 

Our paper is organised as follows. We hrst present the colour decomposition of 
the all-plus, two-loop amplitude, exploiting the multi-peripheral decomposition of the 
underlying tree-level amplitudes. In the next section we describe how the complete 
kinematic structure can be constructed using knowledge from the planar sector and 
tree-level identities. After describing a worked example, we present compact results 
for the full integrand. We perform checks of the universal soft behaviour of the 
amplitude by evaluating the leading 0(e~^) poles of the integrals in the dimensional 
regularisation parameter e. Finally, we draw some conclusions and discuss some 
future directions. 


2 Review of irreducible numerators 


We follow a multi-loop integrand reduction algorithm [51, 53, 54, 58, 72] which uses 
multivariate polynomial division to hnd an integrand representation of the two-loop 
amplitude. This section is intended as a brief overview of the approach; we encourage 
the reader to refer to the literature for more detailed information. 

An integrand-reduced two-loop amplitude has the form 


{p.}) = ia"*" 


Ar({ai}, {pi},-^ 1 ,1*2) 

(2vr)2d 2 ^ A, 4) ’ 


( 2 . 1 ) 


where the sum over runs over graphs T, which are dehned by a specihc set of denom¬ 
inators Da (the set {a} labels the propagators in the graph T.) Associated with each 
graph in the sum is a colour-dressed irreducible numerator Ap; these are functions 
of the external momenta pi, the loop momenta ii and £ 2 , and also of the external 
colour indices ai. Each of these numerators has a colour decomposition 


Ar({ai}, [Pi], £1, £2) = '^(to 

a 


Cr{{ai}) Ar{{pi},£i,£2) , 


( 2 . 2 ) 


where we must explicitly determine the permutation sum a and the associated colour 
factors Cr- We will present an algorithm to hnd a simple colour decomposition for 
our Yang-Mills amplitudes in the next section. 
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An irreducible numerator Ar({pj},^ 2 ) can be written in terms of monomials 
of irreducible scalar products (ISPs). To determine a set of ISPs, we first choose 
a spanning set of momenta to expand the scalar products along the lines of the 
van Neerven-Vermaseren basis [73]. By re-expressing the propagators in terms of 
this spanning set of scalar products, we can see that many can be written as linear 
combinations of propagators and can therefore be removed and pushed down into 
simpler topologies. These scalar products are known as reducible scalar products 
(RSPs.) The remaining scalar products are the ISPs; in general the propagators 
will be quadratic functions of them. To hnd a basis set of ISP monomials, these 
additional quadratic relations are removed using polynomial division with respect to 
a Grobner basis of the relations. This technique has been developed in the public 
code BasisDet [53]. When working in d = 4 — 2e dimensions we also include three 
extra-dimensional ISPs 

= ( 2 - 3 ) 

Once a basis set of ISP monomials is identihed, their rational coefficients are 
computed from the generalized unitarity cuts of the amplitude. As we take all the 
propagators contained in a particular graph on shell, the cut amplitude factorises into 
a product of tree-level amplitudes summed over internal helicity states. Following 
our schematic notation we can write this as 


Cutr = 


■ hi aGF 

r'Dr 


cut(r) 

Ar'({Pi}. 4) Ilaer r*a({Pi}. <i. 4) 


(2.4) 


- cut(r) 


where the trees are those associated with each vertex in the graph P. Making the 
distinction between the cut associated with a graph P and the irreducible numerator 
associated with the same graph is crucial for understanding this construction. The 
irreducible numerator contains only that information which is required on the cut 
associated with P, and which is not captured by irreducible numerators of graphs 
P' that are “larger” than P, in the sense that the propagators contained in P' are a 
proper superset of the propagators contained in P. In other words, by applying the 
cuts in a top-down approach we can isolate each topology systematically subtracting 
the higher-point singularities. 

We will frequently specify the irreducible numerator associated with a graph P 
as A(P) for clarity; one should remember that the function A(P) depends on loop 
and external momenta. Furthermore, throughout this paper we will adopt an index 
notation for the graph labels which lists the number of propagators in each of the 
three two-loop branches £i, G and £i + £ 2 - In addition, we add extra labels to 
distinguish between topologies of this type. We follow the convention that the right 
branch is hrst index, the left branch the second and hnally the central branch in 
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the last entry. For example, the planar pentagon-box reads A 431 = AOiCj), while 
donble box with hve legs in a non-planar ordering is written Asai-s^j = A(M)- 
A complete dictionary between this nomenclatnre and the graphs relevant for onr 
two-loop, hve-point calcnlation is given in Table 1. 

As shown in ref. [58], there is only ever a single branch to the set of solntions to 
the on-shell eqnation in d dimensions, which simplihes the inversion of the system in 
eq. (2.4) to hnd the coefficients of the ISP monomials in Ar, thongh at the cost of 
an increased nnmber of monomials with respect to the fonr-dimensional case. 

There are two important remarks abont this constrnction. The integrand repre¬ 
sentation of eq. ( 2 . 1 ) is not nniqne, and there are different choices for both the set 
of ISPs and the set of monomials. Different sets of spanning vectors will resnlt in 
different ISPs and the polynomial division reqnires a choice of monomial ordering. 

In the following sections we will exploit two important conseqnences of this ap¬ 
proach. Firstly, we will restrict the form of onr irredncible nnmerators to ensnre 
that the choice of ISPs and monomials satisfy the basic symmetries reqnired by onr 
colonr decompositions. Secondly, we will make use of the factorisation of irreducible 
numerators into ordered tree-level amplitudes. These tree amplitudes satisfy a num¬ 
ber of relationships among different orderings. This fact will allow us to determine 
all non-planar cuts of the two-loop all-plus amplitude from the planar irreducible 
numerators computed in ref. [58]. 

3 Colour decomposition 

The main result of this work is the construction of the complete hve-point, two- 
loop, all-plus amplitude in Yang-Mills theory. As we mentioned in Section 2 , it is 
necessary to choose a particular colour decomposition. This decomposition picks a 
set of colour tensors describing the colour structure of the amplitude. At the same 
time, it specihes an associated set of cut diagrams which must be computed. Each of 
these cut diagrams is, in turn, associated with a unique irreducible numerator. Thus 
the colour decomposition that we pick is of central importance, because it determines 
the set of irreducible numerators that we need to calculate. 

In this section, we describe the general algorithm that we used for constructing 
an appropriate colour decomposition of the amplitude, before applying this algorithm 
to the specihc case of the two-loop hve-point amplitude. 

3.1 Multi-peripheral colour decomposition 

Our algorithm is applicable to the general case of an L-loop Yang-Mills amplitude. 
Following the generalized unitarity principle, we begin by writing the amplitude as 
a sum over all colour-dressed cuts. Diagrammatically, these cuts consist of vertices 
formed from colour-dressed tree amplitudes which are joined by on-shell propagators. 
At two loops the set of colour-dressed cuts can be classihed by two basic topologies: 
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a{n — l) 


1 ------ n 

Figure 1: Multi-peripheral diagram for the colour factors in eq. (3.1). 

the genuine two-loop topologies are shown in hgure 2, and the one-loop squared (or 
butterfly) topologies in hgure 3. 

The central idea is to build the loop-level colour decomposition using knowledge 
of the underlying tree-level amplitudes. There are a variety of well-known presenta¬ 
tions of these trees. We hnd it convenient to use the DDM form [70, 71] for the tree 
amplitudes: 

y* CtlCt(T{2)^l J^n-4aa(n-2)^n-3 i>n-3<^(7(n-l)<^n 

cr^Sn — 2 

xA®(l,a(2),... ,cr(n - l),n), (3.1) 

where are proportional to the standard structure constants in SU (W), 

the gauge group of our Yang-Mills theory. The main advantage of this form of the 
amplitude is that it contains (n — 2)! colour structures, as compared to (n — 1)! 
in the standard trace-based decomposition, for example. This fact helps to reduce 
the number of generated diagrams; in particular, an algorithm based directly on the 
trace decomposition of tree amplitudes generates a larger set of diagrams, some of 
which are rather obscure. 

Each of the colour structures in the tree decomposition is a string of group 
theoretic structure constants For an n-gluon amplitude the decomposition is 

constructed by hxing the position of any two gluons at either end of this string. The 
(n — 2)! permutations of the remaining gluons between the ends of the string form 
the set of colour factors each of which is associated with a colour-ordered tree of 
the same ordering. Pictorially, the colour structures look like combs, as shown in 
hgure 1. 

It is straightforward to build the loop colour structure from these DDM tree 
colour structures. The loop structure follows directly from the cut diagram: one 
simply inserts the DDM trees at the vertices; propagator lines connecting trees indi¬ 
cate that the ends of the DDM combs at either end of the propagator have the same 
colour index to be summed over. Notice, however, that we must pick two special lines 
in the DDM form of the tree amplitudes (corresponding to lines 1 and n in hgure 1.) 
These lines are on opposite ends of the DDM colour strings, so one can informally 
think of this choice as picking two lines and “stretching” the colour ordered tree 
between these two ends. We make canonical choices of which legs to pick as special, 
depending on the number of propagators that connect to the three-point amplitude. 
These choices are: 
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Figure 2: Inserting the DDM tree basis into the colour dressed cuts of a two-loop ampli¬ 
tude. The upper insert (a) shows the simple case of two loop propagators, while the lower 
insert (6) shows the case of three loop propagators. The sums run over the permutations 
of the external legs in the tree-level amplitude. 

• Two propagators 

In this case, it is natural to “stretch” using the two propagators as the special 
legs. Thus we build the colour structure by pasting a DDM multi-peripheral 
colour structure between the two propagators. We must sum over every order¬ 
ing of the external legs. Pictorially, the operation is show in the upper insert 
(a) of hgure 2. 

• Three propagators 

In the case of three propagators, we select two out of the three propagators 
to be the special lines in the DDM presentation. Notice that this choice hides 
some of the full symmetry of the diagram. In constructing the DDM tree, 
the propagators we have selected must be at the end of the multi-peripheral 
structure; we must sum over the positions of the other legs. The result is a 
sum of diagrams, as shown in the lower insert (6) of hgure 2. 

• Four propagators 

We again choose two propagators to “stretch” the cut amplitude into a DDM 
tree. At two loops, we only encounter this case in the butterhy topologies. We 
choose upper and lower propagators on the right side of the diagram as special; 
by symmetry, the result is the same as if we chose upper and lower propagators 
on the left of the diagram. The insert of hgure 3 sketches out the procedure. 

In this way, we build a set of colour structures. The kinematic structure associ¬ 
ated with each colour structure is easily understood. Each time we insert a particular 
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Figure 3; Inserting the DDM tree basis into a colour dressed cut of a butterfly topology 
at two-loops. There are four loop propagators in this case, and the insert (a) shows the 
result of inserting the DDM tree decomposition fixing the two right legs. The sums run 
over the permutations of the external legs in the tree-level amplitude. 


DDM colour trace, we also pick up a factor of the associated colour-ordered tree am¬ 
plitude. Thus, the orientation of the legs in the kinematic diagram, associated to an 
irreducible numerator, is the same as in the colour structure; of course, the “stretch¬ 
ing” procedure does not produce new propagators in the irreducible numerator. 

One advantage of using the DDM presentation of the amplitude at tree level 
is that the KK relations are automatically satished. Our procedure recycles this 
property to loop level: we automatically generate a set of colour diagrams that is 
KK-independent. Along the way, we generate ordered diagrams for the kinematics. 
The same procedure works at L loops; the amplitude is expressed as 



^ jfl 

KK-independent j=l 
IPI graphs Ti 


d% lQA,{i) 
{2nYS, D,{i) 


(3.2) 


where Si are the symmetry factors of the graphs, Di denote the products of the 
(inverse scalar) propagators, and the Aj are irreducible numerators for appropriate 
colour factors Cj generated through our algorithm. Now let us see this procedure at 
work in the context of the hve-point, two-loop amplitude, which is our main focus. 


3.2 Five-point, two-loop amplitude 

Now we describe the colour structure of the hve-point, two-loop amplitude. We 
concentrate on the diagrams that do not vanish in the all-plus case according to 
ref. [58] and our calculations in section 4. A generic hve-point two-loop amplitude 
can be constructed by a straightforward extension of the present discussion. 

Let us hrst write the amplitude and then explain its content in more detail. We 
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label Aj and their colour factors by their diagrams directly in the formula: 


^(2)(i,2,3,4,5) = 


o-eSs 


c 




+ A 


- 


+c 




m: 

+ ... 


A- 


4 


where the integration operator I acts on every Aj as 


I [A.] ^ 


d‘‘ljd‘‘l2 Aj 
(2ir)“ A' 


(3.4) 


The explicit symmetry factors compensate for the over-counts introduced by the 
overall sum over permutations of external legs. For convenience, we recapitulate all 
these Aj in table 1 , where for each irreducible numerator Aj we also show its diagram, 
that of its colour factor, as well as the set of its non-equivalent permutations. 

The first three graphs in table 1, A 431 , A 332 and A 422 , are the master diagrams 
corresponding to the maximal cuts. They are purely trivalent, thus their colour 
factors are unambiguously dehned by their proper graphs. 

The next three graphs, A 331 ;A 4 -j, A 232 ;Mi and A 322 ;Mi, have a four-point vertex 
with two external and two internal edges. The two external legs automatically enter 
in the permutation sum with two possible orderings, hence multi-peripheral sub¬ 
graphs are naturally obtained by hxing the internal lines, as in the insert (a) of 
figure 2. “Stretching” the four-point vertex by these lines gives a master graph for 
each colour factor. 

The following two diagrams, A 33 i. 5 ij and A 33 i. 52 , 2 , share the same graph struc¬ 
ture, up to the ordering of the four-point vertex. The apparent asymmetry introduced 
by our selecting these two diagrams, and omitting the graph with the external leg 
in the right loop, is an artefact of our colour decomposition. One could make other 
choices; the KK relations satished by the trees and their symmetries ensure that any 
choice is valid. 

To expand the four-point vertex in the ninth graph, A 322 ; 5 Li , we hxed the internal 
lines of the “diamond” subdiagram, hence its colour factor is 0 * 332 , but with a minus 
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Numerators 

Graphs 

Colour factors 

Permutation sum 

A 431 



S's/Vertical flip 

A 332 

\ y 

y k 

Sj__ y 

y —— k 

S' 5 /Vertical & horizontal flip 

A 422 



S's/Vertical & diamond flip 

A331;Mi 

:=A 


^5 

A232;Mi 


N,__ y 

y —— k 

S's/Vertical flip 

A322;Mi 



S's/Diamond flip 

A331;5Li 

y — 1 —k 


S's/Horizontal flip 

A331;5L2 


N|__ y 

y —— k 

S's/Horizontal flip 

A322;5Li 



^5 

A 43 O 



S's/Vertical flip 

A330;Mi 

M 


^5 

A330;5Li 

M 


S's/Horizontal flip 

A330;5L2 

M 

S|__ y 

y —— k 

S's/Horizontal & vertical flip 


Table 1: The irreducible numerators that are nonzero for the all-plus five-point two-loop 
amplitude, along with their colour factors and reduced permutation sums. 


sign due to one flipped vertex. The other permutation of the four-point vertex 
corresponds to the same topology and is present in the overall permutation sum 
with the right permuted colour diagram. 

The colour factor of the planar graph A 430 follows in a straightforward manner 
from our algorithm (see hgure 3), yielding C 431 as its colour factor. Its descendant 
A 330 ;Mi is more interesting, since it is the only graph in the all-plus case with two 
four-point vertices. They can be treated independently by linearity of colour decom- 
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position. The external one is thus expanded in the same way as in A 331 .M 1 , producing 
C 430 as an intermediate step. Expanding the internal four-point vertex gives C 431 
again. 

The last two graphs, A 330 ; 5 Li and share a hve-point vertex. To explain 

their colour factors, let us consider the corresponding colour-dressed cut: 


Cut 


330;5L — 


= c 


m: 


c 


Cut 




><;) ( 3 . 5 ) 

m:)' 


We obtain the multi-peripheral decomposition of the hve-point vertex by hxing the 
two right-hand loop edges and permuting the other three edges. The graphs in the 
second line can be vertically hipped to put leg 3 downstairs to match the presentation 
in table 1. Obviously, an equivalent decomposition could be achieved by hxing the 
loop edges on the left, which would change the orientation of leg 3 in the superhcially 
non-planar graphs A 330 ; 5 i, 2 - The S' 5 -summation in eq. (3.3) ehectively symmetrises 
the colour structure over the two choices of multi-peripheral decompositions. 

In the present paper we can avoid lower levels of the graph hierarchy thanks to the 
simplicity of the fully symmetric helicity conhguration, but it already incorporates 
all the key elements of the general colour structure. 


4 Kinematic structure 

With our colour decomposition at hand, we turn our attention to the kinematic 
structure of the amplitude. We need to compute an irreducible numerator associated 
to each diagram in eq. (3.3); as Frellesvig, Zhang and one of the current authors 
have already computed all the planar irreducible numerators [58], our task is to 
determine the remaining non-planar numerators. Of course, these numerators can 
be computed directly from their cuts. However, as we will see, it is easy to determine 
the complete set of non-planar irreducible numerators for this amplitude from its 
planar numerators and the knowledge of tree-level amplitude relations. 

4.1 Non-planar from planar 

The non-planar numerator A 332 = A()eiii() can, of course, be obtained directly from 
its cut. However, we can avoid calculating this non-planar cut explicitly by relating 
it to a planar cut. We do so in two steps: hrst, we coalesce two (ordered) three-point 
trees into a limit of an ordered four-point tree; then we use the BCJ relations [59] 
satished by the ordered four-point tree to reorder the legs until the complete diagram 
becomes planar. 
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In more detail, we use the following well-known relation, which is satisfied by 
on-shell amplitudes in the cuts: 


kl(°)(l,2,-(l + 2))kl(°)(l + 2,3,4) = {si3 2l(°)(l,3,2,4)}| 


Si2=0’ 


4 1 


= < Sl3 


X 


(4.1) 


Sl2=0 


where Sij = (pi + Pj^ are the standard Mandelstam invariants. Since this identity 
is of central importance for us, we present a short proof. A four-point tree ampli¬ 
tude can be constructed from two three-point amplitudes using the BCFW recursion 
relation [74, 75]: 


2, 3,4) = —2, -(l+2))X°)(l + 2,3,4), (4.2) 

Sl2 


where hat signs on the right-hand side indicate that the amplitudes are evaluated on 
complex kinematics for some BCFW shift of external legs. The exact complex value 
of the shifted internal momentum ( 1 - 1 - 2 ) is dehned by the on-shell condition 


Si2 — Si2 + 0!Z — 0, 


(4.3) 


where the precise expression for a depends on the particular BCFW shift. The key 
point is that Si 2 is a linear function of z, with the property that in the limit Si 2 —?• 0 , 
z —?• 0. In this limit eq. (4.2) becomes 


{si2A(°)(1,2,3,4)} 


A(°)(1, 2, -(1 + 2))X°)(1+2, 3,4). 


(4.4) 


Notice that the left-hand side contains a nonzero contribution due to the pole in Si 2 . 
Now we can remove the factor of Si 2 on the left-hand side of eq. (4.4) by using the 
BCJ amplitude relation [59], 


si2A(°)(1,2,3,4) = si3A(°)(1,3,2,4). (4.5) 


This proves the identity (4.1). 

We proceed by applying our identity (4.1) to tree amplitudes inside the non- 
planar cut, rearranging the diagram until it becomes planar. It is simplest to begin 
with maximal diagrams, and then to continue to topologies with fewer propagators. 
We will work through the calculation of A 332 , displayed in Figure 4, as an example; 
we computed all non-planar irreducible numerators using the same technique. 
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5 


1 


^2 



4 2 

(a) A332(123 45,4,^2) (b) A422(123 45,£ 2 ) 

Figure 4; The two non-planar maximal topologies. 
The calculation starts at the level of the cuts: 


Cut332 — 



where = Pi + - ■ ■ -hpj, and we understand that all internal helicities are implicitly 
summed over while all exposed propagators are cut. These cuts are decomposed into 
irreducible numerators as 


Cut332 — A332(123 45 ,£i,£2), 

Cut33i;5Lj = A33i;5i;,j (12345, £ 1 , £ 2 ) + 


A43i(12345,£i,£2) 


(^1 -T123)^ 

A 431 (34512, P 345 — i2,Pi2 — ^ 1 ) 


(4.7a) 

(4.7b) 


(4 - P345)^ 

Using the fact that {i\ — ^ 123 )^ = —{^2 — 7 * 345 )^ on this cut, we see that 

A332(12345, £ 1 , £ 2 ) = (^1 — P123)^A 331;5Li (12345,A,4) + A43 i(12345,£i,£2) 

~A 43 i (34512, P 345 — f' 2 ,Ti 2 — f*i)- 


(4.8) 


A similar calculation for A 422 leads to 


A 422 (12345, A, 4) = A43i(12345, 7i, 4). (4.9) 

So far the obtained non-planar numerators are valid only on their cuts, but they 
can be extended off shell. We may simply express the numerators in terms of a 
given set of ISPs and then dehne off-shell numerators unambiguously through these 
ISPs. The value of a given numerator depends on the choice of ISP basis off shell (in 
contrast to the situation on shell, of course). In this way, we determine a valid set 
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Graphs 

ISPs 

RSPs 

A332 

4 ■ (P5 -P 4 ), 

(^1 -Pl)^ (4 -Pl2)^ 


^2 ■ (Pl -P2), 

^1, (^2 -P5)^ (4 -P45)^ 


{il - £2) ■ P 3 

(G+ 4 +P 3 )^ (G+ 4 )' 

A422 

(£1 + 2^2) • (Pl -P3), 

(^1 -Pl)^ (G -Pl2)^ (G -Pl 23 )^ 


(£i + 2 £ 2 )-P 2 , 

to to 

to 

1 

o\ 

to 


h ■ (P 5 - P 4 ) 

{£1 + £2+ Pi)'^, {£1 + £2)"^ 


Table 2; The choices of ISPs and RSPs for the two non-planar masters, where the RSPs 
are chosen as the propagators of the respective graphs. Additionally, the higher-dimensional 
ISPs fjLij are shared by all topologies. 

of non-planar irreducible numerators. Notice that the ISP monomial choices made 
in the planar sector, such as the higher powers of /i^ preferred over high powers of 
(ii-pj), will then be easily translated to the non-planar numerators. 

It is very useful to maintain the symmetries of the underlying graphs in this 
off-shell continuation. We achieve this by choosing an appropriate basis of ISPs on 
a graph-by-graph basis. One engineers the ISP basis such that the loop momentum- 
dependence of each irreducible numerator is captured by a set of ISPs, which map 
into one another under the graph symmetries without using any cut conditions. The 
symmetries of the maximal non-planar graphs, for example, are 

A332(12345,4,4) = -A332(21354,pi2 - 0 ,P 45 -4) = A332(54 321,0,0), (4.10a) 

A422(123 45, ii, £ 2 ) = —A 422 (32145, P 123 — ~ 0) = A422(123 54, ii, —ii — O)- 

(4.10b) 

These symmetries motivate our choices of ISPs, given in Table 2. For instance, 
the second symmetry of A 332 in (4.10a) leads to a map of ISPs 

h ■ (Ps - Pi) ^ 4 ■ (Pi -P 2 ), 

(£ 1 -4)-Ps ^-(4-4)-Ps- (4.11) 

After we express loop-momentum dependence in eqs. (4.8) and (4.9) in terms of 
the ISPs of Table 2, using the fact that the RSPs (cut propagators) are zero on 
shell, we are left with appropriate off-shell irreducible numerators. These are listed 
in Section 5. Note that the function A 422 is not the same as the function A 431 
despite the on-shell equation (4.9): different ISPs are chosen to make different off- 
shell symmetries manifest. 

We obtained irreducible numerators for lower-level non-planar diagrams in the 
same way, using the BCJ relations on cuts and extending the results off shell. For 
the all-plus amplitude at hand we hnd that many lower-level irreducible numera¬ 
tors vanish. In other words, the higher-level numerators capture the unitarity cut 
structure of the full amplitude, which is given below. 
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5 The full-colour five-gluon all-plus amplitude 

In this section we present a complete snmmary of all kinematic nnmerators con- 
tribnting to the colour decomposition in eq. (3.3). We include both planar [58] and 
non-planar irreducible numerators computed using the technique described in Sec¬ 
tion 4.1. The result is presented unrenormalized including the dependence on the 
spin dimension Dg of the gluon, which is equal to 4 in the FDH scheme and 4 — 2e 
in CDR [37]. The dependence on the extra dimensional ISPs fiij = i\ • Sj can 
be collected into three general functions, 

Fi = {Dg —2) infill ^22 + (hii + 7 ^ 22 )^ + 2(/rii -|- /i 22 )/^i 2 ) + 16(/i^2 ~ ^ 11 /^ 22 ), (5.1a) 
F 2 = 4(Zls —2)(/iii-|-/i 22 )/^i 2 , (5.1b) 

F 3 = {Dg — 2 )^fj,iifi 22 - (5.1c) 


The remaining coefficients are expressed using the standard spinor-helicity formalism. 
In particular, we denote 


trs = ‘iie^^pap’(p2P3Pl 

= tr+(1234) - tr_(1234) (5.2) 

= [12] (23) (34) (41)-(12) [23] (34) [41], 


since tr± = tr (4(1 ± 75 ) 1234 ). We also make use of “spurious” directions i 
to find compact representations of the integrands. 


m 


^abc = 


{bc)[ca] {ac)[cb] (6|7^|a] 


^ab 


^ab 


order 


(5.3) 


The full amplitude reads 

^(2)(i+,2+,3+,4+,5+) ^ 




c 







iA(;>^)+A(;)xr“)+iA(w; 


+c 



(5.4) 


4^yA><: 




4 - 
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The three maximal graphs are 


A4.cii = A( j = 


'Sl2'S23'S45-^l 


A 332 — A 


5 ^ -^2 ^1^1 


(12)(23)(34)(45)(51)tr5 

'Sl2'S45-^l 


(tr_|_(1345)(£i + + 'Si5'S34’S45) ) 


A 422 — A 




4(12)(23)(34)(45)(51)tr5 

X (^S23tr+(1345)(2si2 - 4£i-(p5 -^ 4 ) + 2(A - ^2)-Pz) 

- S34tr+(1235)(2s45 -4£2-(pi - P2) - 2(£i - i2)-P3) 

— 4s23<S34Si5(£i — £ 2) 'Ps ^ ) 

’Si2'S23'S45-^1 


(12)(23)(34)(45)(51)tr5 

X ^tr+(1345)^^l-(P5 - Pi) - + Xi5S34S45^. 


(5.5) 


Meanwhile, the graphs at level 1 are 


'5^2 *1| 

A43n = A ( 


Si2tr+(1345) 


(2(£i-a;i23) + "S 23 ) 


‘^331;5Li — ^ 


'^331;5L2 “ ^ ( 13 


^322;5Li — A 


A331;Mi — A 




2(12)(23){34>{45>(51)s, 3 
X A + F3A±h)!±il5 

\ >545 

>S12>S23’S34>S45'S51.Fi 

(12)(23)(34)(45)(51)tr5’ 

_ >S12>S45-^1 

“4(12)(23)(34)(45)(51)tr5 
X ( 2 S 23 -S 34 -S 15 - S23tr+(1345) + S34tr+(1235)), 


t 


2(12)(23)(34)(45)(51)tr5 
X (s23S45tr+(1435) - si5S34tr+(2453)), 

= A322;Mi = A ^ = A232;Mi ~ 

S34S45tr+(1235)Fi 


(12)(23)(34)(45){51)tr5’ 


(5.6) 
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Finally, the graphs at level 2 are 


A 


330; Ml 


= A 


^330;5Li — ^ 
1 


:r: 


(g45 - gi2)tr+(1345) 
■2(12)(23)(34)(45)(51 )si3 
1 


F 2 + F; 


{^l + £ 2 )^ + -545 


"545 


X <1 - ( tr+(1245) - 


X F 2 + Fa 


(12)(23)(34)(45)(51) 
tr+(1345)tr+(1235)\ 

'Sl3'S35 / 

4(^1-ps) (^2 '^ 3 ) + (^1 + ^2)^("Si2 + < 545 ) + Sl2S45\ 
'Sl2'S45 / 


(£1 + ^' 2)^515 (5.7) 

+ tr+(1235) { AlA! -'^I2l(l + _p,)2') 

V 2SqK Si 9 V Sir S.9RS/1R / 


tr+(1345) 


\ 2s35 S12 V S35 <S35S45 

(^1+^2)^ ^ 2 ' Pz (^ , 2(£i-a;i23) , >545 — S12 


2 s 


13 


A 


330;5L2 


= A 


X 


5 V ^2 ^1 y 1 

M, 


>545 

7^3 


1 + 


>513 


+ 


>5i2>5i3 


(7i -Pi)' 


2(12)(23)(34)(45)(51)si2 


(S 45 - si2)tr+(1245) - ('tr+(1245) - tr+(1345)tr+(1235) \ 

\ >513>535 / 


S45tr+(1235) 

>535 

Si2tr+(1345) 

>513 


2(t'2-lX’543) H- ^(£2 — Pb)‘ 

>545 

2(t'i-x;i23) H —— -—(7i — Pi)" 

>512 


We have found a representation of the full amplitude with no topologies with 
fewer than six propagators. We note that there are nonzero cuts at the integrand 
level, but the resulting integrals are scaleless and hence zero in dimensional regular- 
isation. We have checked additional cuts at levels 2 and 3 to ensure that no nonzero 
topologies remain. To find an integrand with this property, the ISPs (f'i-a;i 23 ) and 
(£2 ■ 1 ^ 543 ) in the numerators A 330 ; 5 L 2 and A 330 ; 5 Li were upgraded to include terms 
proportional to (£1 — Pi)^ and (£2 — P 5 )^> 


6 Checking the soft divergences 

Since the all-plus helicity conhguration is zero at tree level, the universal infrared 
(IR) structure is the same as that of an ordinary one-loop amplitude. The poles of 
our two-loop amplitude should therefore be equivalent to those of the hnite one-loop 
amplitude multiplied by the infrared pole operator including the sum over colour 
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correlations [76] 


^(2)(i+2+3+4+5+) = 

Z ri-T,o^W(l+ 2 +, 3 +, 4 +, 5 +) + 0 (£-'), 


( 6 . 1 ) 


j¥=i 


where we have used the o symbol to indicate that the colour matrices (in this purely 
gluonic case they will all be structure constants) should be inserted into the colour 
factors of the one-loop amplitude. The standard loop prefactor is given by ^ 


r(i + £)r^(i-6) 

(47r)2-<r(l-2e)' 


( 6 . 2 ) 


The one-loop amplitude to all orders in e can be found in ref. [68]. 

There are a number of difficulties in checking eq. (6.1) in full since the five-point 
planar and non-planar integrals required are still unknown at this time. Resorting 
to numerical evaluation, as has been done in the planar case [58], is computationally 
prohibitive for two reasons. Firstly, the full colour expansion contains a large num¬ 
ber of dimension-shifted integrals (~ 0(1000)) - an order of magnitude more than 
the leading colour terms. Secondly, there is no Euclidean region for the complete 
amplitude, and so contour deformation must be performed for many of these inte¬ 
grals, making them more complicated than the planar cases. This task is probably 
achievable with public tools like FIESTA [77] and SecDec [78]. 

There is, however, a much simpler method to check the leading soft singularities 
up to which can be done analytically. In the leading soft limit, the colour 

correlations drop out of eq. (6.1): 


^(2) (1+2+3+4+5+) = 2+,3+,4+,5+) + C>(e"^). (6.3) 

Clearly this is a weaker check than the full IR poles, but it does require non-trivial 
properties of the non-planar sector. The butterfly (one-loop squared) topologies, are 
all finite and therefore not relevant for the IR properties. Scattering amplitudes in 
the soft or eikonal limit have many remarkable structures and universal properties. 
The interested reader may like to turn to ref. [79] for a recent introduction to the 
subject. 


6.1 Evaluating the massless double box in the soft limit 

The extra simplicity in our all-plus loop amplitude that sets it aside from most 
two-loop amplitudes is that the integrals contain at most a single soft divergence, 

^The normalisation of the integrals in this section is different by a factor of per inte¬ 

gration with respect to the default choices in FIESTA and SecDec. 
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rather than the maximum double soft divergence. We can therefore break our loop 
amplitudes up into sums of regions with soft singularities and evaluate the amplitude 
in the limit. In this limit the integral factorises into a product of two one-loop 
integrals and can be evaluated to extract the leading 0{e~^) divergence. 

All of the poles of our amplitude are contained in the topologies proportional to 
the same dimension shifting numerator: 

Fi = {Dg — 2)(3p.iip.22 + hii + f ^22 + 2/ii2(/iii -|- p- 22 )) + 16(/r^2 ~ f^nfJ‘ 22 ), (6-4) 

which also has a simple behaviour, 

limFi = {Dg — 2)/i22- (6.5) 

£1-5-0 

Taking the example of the two-loop double box, we find two soft regions by 
taking the limit of either loop. We find a soft singularity whenever we have two 
adjacent massless legs in one of the loop integrations. In each case, we factorise into 
an IR divergent triangle and a dimension-shifted box: 


w ^ :k: i-y k: • 


(6.6a) 


('■)□□('') in] (D. - 2)/*-=- r)>=) /*-=■ (‘)l( j Ia?,]. (6.6b) 


Recalling the one-loop integrals. 


r4-2e 




^n] 


(dvr) 2 6 




' ' (47r)2 Si2e2 

and summing the two regions, we quickly arrive at the result: 


+ 0(£-), 


(6.7) 

( 6 . 8 ) 


r4-2e 




1 Ds-2 1 , 

+ 0{e-^). 


(47r)^ 35126^ 


(6.9) 


6.2 Soft divergences of the five-point integrals 

By following the method described in the previous section, we have derived the com¬ 
plete set of integrals required for the 0(€~^) part of the amplitude. All of the integrals 
have been checked numerically using the sector decomposition methods implemented 
in FIESTA [77] and SecDec [78]. Some of these integrals have been computed long 
ago in 4 — 2e dimensions and can be used to write the full integrals including finite 
terms via the dimensional reduction identities implemented in LiteRed [80] and 
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IBP relations from FIRES [81].^ We have performed this task for the planar donble 
box with an off-shell leg. Thus we arrived at the following soft limits for the integrals: 


r4-2e 


:^ 3 [ni 


1 D.-2 






(dvr)^ 3si2S23e^ 

1 — 2)(2si5-|-S 25 ) 


(dvr) 


12si2S23e^ 


(6.10a) 


+ C>(e-'), (6.10b) 


j4-2e 

j4-2e 

j4-2e 


4-^ ^2 

5v -^2 


4^ ^2 

5^ -^2 


[ni = o(e-‘), 

[n (fr(P5 -P4))l = 0(e"‘), 

[n(((’.-«-P3)l = 0(e-‘), 


r4-2e 




1 D ,-2 


+ Oie-% 


(6.11a) 

(6.11b) 

(6.11c) 

(6.12a) 


(d7r)4 3si2S23e^ 

-P4))l = + ^34-333) ^ 

(6.12b) 




(dvr) 


125125236 ^ 


r4-2e 


[U] = 74^1%^ (- + -)+ (6.13) 


(d 7 r )4 6 e 2 \si2 S45 


r4-2e 


r4-2e 


j4-2.^ 


^;)in] 


1 Ds-2 

(dvr)^ 65126^ 






1 - 


1 Ds-^ 

(dvr)"^ 65456 ^ 




(6.1d) 


(6.15) 


(6.16) 


(6.17) 


Using these results, we have checked that eq. (6.3) does hold as expected for our 
amplitude (S.d). 


)hx;) mi 


1 Ds -2 
(dvr)^ 65456^ 


+ 0(4-‘), 




^We thank Claude Duhr for providing his own computation of the integrals for e"'"e —> 3j 

[82, 83]. 
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7 Conclusions 


In this paper we have explored the impact of tree-level amplitude relations in multi¬ 
loop integrand computations. There were two major aspects to our work. Firstly, we 
exploited the Kleiss-Kuijf relations to find a compact colour decomposition for the 
two-loop amplitude in terms of multi-peripheral colour factors in an analogous way to 
the tree-level and one-loop decompositions of Del Duca, Dixon and Maltoni [70, 71]. 

Secondly, we applied the BCJ amplitude relations [59] to relate all non-planar 
generalized unitarity cuts to the previously computed planar ones. This allowed us to 
easily generate a compact representation the full colour two-loop, five-gluon, all-plus 
integrand building on previous planar work [58]. The soft infrared poles of the full 
amplitude were checked against the well-known universal pole structure. 

We hope that the computational methods developed here will be of good use 
in the necessary extension to more general helicity conhgurations and other 2 —)■ 3 
scattering processes at two loops. They highlight some advantages of relating two- 
loop integrands to tree-level amplitudes via generalized unitarity cuts. As well as 
avoiding the large intermediate steps that make Feynman diagram computations at 
this loop order and multiplicity extremely computationally intensive, we are able to 
build known on-shell symmetries and relations into the amplitude by construction. 

Another interesting aspect of the all-plus amplitude is the continuing connection 
to the previously known amplitudes in A/^ = 4 sYM. Though the dimension shifting 
relation observed at one loop no longer holds, the integrands of our full all-plus 
amplitude and the expressions of Carrasco and Johansson [63] are related by the 
same dimension shifting operator seen in the planar case. For example, we hnd that 


A 


2+,3+,4+,5+) = ^22, ^il2) 2- 3+ 4+ 5+)^ ^ ^ q . (7.1) 


The one-loop squared topologies have the form 


AxyO;T(l^, 2"'", 3^, d"*", 5"'') — (F"i (/in,/i22,/^12) — (/in,/i22, —/X 12 )) A({pj}, £ 1 , £ 2 ) 

+ -^3(7^11) 1 ^ 22 , hl 2 )-B({Pi}, il, -^2); 


(7.2) 


where A and B are some functions of the external kinematics and loop momenta. 
The second term is proportional to {Dg — 2)^ and is a genuine contribution in QCD 
not related to A/" = 4. This additional numerator structure is enough to make the 
off-shell BCJ symmetries non-trivial to satisfy, even though the JV = 4 integrand has 
been cast in such a form. It is an interesting question as to whether this would be 
possible for the amplitudes presented here and one that we intend to explore in the 
future. 
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